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Abstract 



We formulate an effective Schrodinger wave equation describing the quan- 
tum dynamics of a system of DO-branes by applying the Wilson renormal- 
ization group equation to the worldsheet partition function of a deformed 
<7-model describing the system, which includes the quantum recoil due to the 
exchange of string states between the individual D-particles. We arrive at 
an effective Fokker-Planck equation for the probability density with diffusion 
coefficient determined by the total kinetic energy of the recoiling system. We 
use Galilean invariance of the system to show that there are three possible 
solutions of the associated non-linear Schrodinger equation depending on the 
strength of the open string interactions among the D-particles. When the open 
string energies are small compared to the total kinetic energy of the system, 
the solutions are governed by freely-propagating solitary waves. When the 
string coupling constant reaches a dynamically determined critical value, the 
system is described by minimal uncertainty wavepackets which describe the 
smearing of the D-particle coordinates due to the distortion of the surrounding 
spacetime from the string interactions. For strong string interactions, bound 
state solutions exist with effective mass determined by an energy-dependent 
shift of the static BPS mass of the DO-branes. 
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I. INTRODUCTION AND SUMMARY 



Despite the enormous amount of activity over the past few years towards understanding 
the dynamics of Dirichlet p-branes Q, the problem of demonstrating that a system of iV 
moving D-particles can form a bound state is still unresolved. It is relevant to the Matrix 
Theory conjecture in which these particles are interpreted as Kaluza-Klein modes of 11- 
dimensional M-Theory compactified on a circle and are described by the supersymmetric 
N x N matrix quantum mechanics that is obtained from dimensional reduction of M = 1 
supersymmetric Yang-Mills theory in ten dimensions. The existence of such a tower of states 
is equivalent to the statement that the quantum mechanics admits exactly one bound state 
for each N. The original threshold bound state problem was addressed in and recent 
progress has been made in M. However, beyond the proof of existence in the case N = 2, the 
general case N > 2 remains in large part an open problem. In this paper we shall address 
the bound state problem for a system of N non-relativistic, recoiling DO-branes by studying 
their moduli space dynamics in certain limits. 

The effective worldvolume dynamics of a single Dp-brane coupled to a worldvolume gauge 
field and to background supergravity fields is described by the action Jl| 

Id p = %J d p+1 a e-+ y j-^[G aP + B aP + 2'KQ!F a p] +T P J d p+1 a 

(1.1) 

The first term in ( |1.1|) is the Dirac-Born-Infeld action with T p the p-brane tension, a' the 
string Regge slope, <fi the dilaton field, F = dA the worldvolume field strength, and G and B 
the pull-backs of the target space metric and Neveu-Schwarz two-form fields, respectively, to 
the Dp-brane worldvolume. It is a generalization of the geometric volume of the brane tra- 
jectory. The second term is the Wess-Zumino action (restricted to its p + 1-form component) 
with C the pullback of the sum over all electric and magnetic Ramond-Ramond (RR) form 
potentials and Q a geometrical factor accounting for the possible non-trivial curvature of the 
tangent and normal bundles to the p-brane worldvolume. It describes the coupling of the 
Dp-brane to the supergravity RR p + 1-form fields as well as to the topological charge of the 
worldvolume gauge field and to the worldvolume gravitational connections. The fermionic 
completion of the action fll.lj ), compatible with spacetime supersymmetry and worldvolume 
^-symmetry, has been described in ||. For a recent review of the Born-Infeld action and its 
various extensions in superstring theory, see ||. 

While the generalization of the Wess-Zumino Lagrangian to multiple Dp-branes is obvious 
(one simply traces over the worldvolume gauge group in the fundamental representation), 
the complete form of the non-abelian Born-Infeld action is not known. In it was proposed 
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that the background independent terms can be derived using T-duality from a 9-brane action 
obtained from the corresponding abelian version by symmetrizing all gauge group traces in 
the vector representation ||. A direct calculation of the leading terms in a weak supergravity 
background has been calculated using Matrix Theory methods in ||. Based on the Type I 
formulation, i.e. by viewing a D-particle in the Neumann picture and imposing T-duality 
as a functional canonical transformation in the string path integral |J , the effective moduli 
space Lagrangian was derived in [IU| and shown to coincide (to leading orders in a velocity 



expansion) with the non-abelian Born-Infeld action of 0. In the following we will use this 
moduli space approach to D-brane dynamics to describe some properties of the multiple 
D-brane wavefunction. 

The novel aspect of the approach of is that the moduli space dynamics induces 



an effective target space geometry for the D-branes which contains information about the 
short-distance spacetime structure probed by multiple D-particles. Based on this feature, 
string-modified spacetime and phase space uncertainty relations can be derived and thereby 
represent a proper quantization of the noncommutative spacetime seen by low-energy D- 



particle probes [TT[. The crucial property of the derivation is the incorporation of proper 
recoil operators for the D-branes and the short open string excitations connecting them. The 
smearing of the spacetime coordinates y? (in general % = 1, ... ,9 — p label the transverse 
coordinates of the Dp-brane and a — 1, . . . , N the component branes of the multiple D-brane 
configuration) of a given D-particle as a result from its open string interactions with other 
branes can be seen directly from the formula for the variance 

(Ay?) 2 = [(Yi - Yri N f] aa = £ \ Y f\ 2 (1-2) 

b^a 

where Yf b are the w(iV)-valued positions of the D-particles (a = b) and of the open strings 
connecting branes a and b (a ^ b), and In is the N xN identity matrix. The recoil operators 
give a relevant deformation of the conformal field theory describing free open strings, and 
thus lead to non-trivial renormalization group flows on the moduli space of coupling con- 
stants. The moduli space dynamics is thereby governed by the Zamolodchikov metric and 
the associated C-theorem. Physically, the recoil operators describe the appropriate change of 
quantum state of the D-brane background after the emission or absorption of open or closed 
strings. They are a necessary ingredient in the description of multiple D-brane dynamics, 
in which coincident branes interact with each other via the exchange of open string states. 



The quantum uncertainties derived in ]10|jn[ were found to exhibit quantum decoherence 



effects through their dependence on the recoil energies of the system of D-particles. This 
suggests that the appropriate quantum dynamics of DO-branes should be described by some 
sort of stochastic string field theory involving a Fokker-Planck Hamiltonian. 
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As in H[7|], the derivation in |10| assumes constant background supergravity fields. How- 
ever, another important ingredient missing in the moduli space description are the appro- 
priate residual fermionic terms from the supersymmetry of the initial static D-brane config- 
uration. While the recoil of the D-branes breaks supersymmetry, it is necessary to include 
these terms to have a complete description of the stability of the D-particle bound state. 
As shown in JL2|j , the energy of the bound states of D-branes and strings is determined by 
the central charge of the corresponding spacetime supersymmetry algebra. Nonetheless, the 
bosonic formalism that we display below can be exploited to a large extent to describe at 
least heuristically the quantum phase structure of the multiple D-particle system and, in 
particular, determine the mass and stability conditions of the candidate bound state. One 
reason that this approach is expected to yield reliable results is that we view the system 
of D-branes and strings as a quantum mechanical system (rather than a quantum field 
theoretical system as might be the case from the fact that T-duality is used to effectively 
integrate over the transverse coordinates of the branes), with the D-brane recoil constituting 
an excitation of this system. The recoiling system of D-branes and strings can be viewed 
as an excited state of a supersymmetric (static) vacuum configuration. The breaking of 
target space supersymmetry by the excited state of the system may thereby constitute a 
symmetry obstruction situation in the spirit of [[U| . According to the symmetry obstruction 
hypothesis, the ground state of a system of (static) strings and D-branes is a BPS state, but 
the excited (recoiling) states do not respect the supersymmetry due to quantum diffusion 
and other effects. Phenomenologically, the supersymmetry breaking induced by the excited 
system of recoiling D-particles will distort the spacetime surrounding them and may result 
in a decohering spacetime foam, on which low energy (point-like) excitations live. This mo- 
tivates the study of non-supersymmetric D-branes recoiling under the exchange of strings. 
Such quantum mechanical systems exhibit diffusion and may be viewed as non-equilibrium 
(open) quantum systems, with the non-equilibrium state being related naturally to the pic- 
ture of viewing the recoiling D-brane system as an excited state of some (non-perturbative) 
supersymmetric D-brane vacuum configuration. 

The main relationship we shall exploit in obtaining the quantum dynamics of multiple 
D-particle systems is that between the Dirichlet partition function in the background of 
Type II string fields and the semi-classical (Euclidean) wavefunctional ty\Y l ] of a Dp-brane. 



This relation is usually expressed as [14, 



Z = I DY l m\Y l ] (1.3) 

The wavefunction ty[Y l ] is expressed in terms of the generating functional which sums up all 
one-particle irreducible connected worldsheet diagrams whose boundaries are mapped onto 
the D-brane worldvolume. Integration over the worldvolume gauge field is implicit in \1/ to 
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ensure Type II winding number conservation. Dirichlet string perturbation theory yields 

oo 

nY i } = ex P J2^ h ~ 2) ^h[Y l ] (1.4) 



h=l 



where Sh denotes the amplitude with h holes, in which an implicit sum over handles is as- 
sumed. However, as we will discuss in the following, the identification ( |1.3|) is not the only 
one consistent with the approach to D-brane dynamics advocated in ||10|| , and one may in- 
stead identify the worldsheet Dirichlet partition function, summed over all genera, with the 
probability distribution corresponding to the wavefunction Using this identification, the 
Wilson renormalization group equation has been proposed as a defining principle for obtain- 
ing string field equations of motion, including the appropriate Fischler-Susskind mechanism 



for the contributions from higher genera |TJ[]. When applied to Dirichlet string theory, we 
shall find that the consistent D-brane equation of motion follows from the renormalization 
group equation. 

More precisely, within the framework of a perturbative logarithmic conformal field the- 
ory approach to multiple D-brane dynamics ||10|| , we will show that the intricate quantum 
dynamics of a system of interacting D-particles is described by a non-linear Schrodinger 
wave equation. The corresponding probability density is of the Fokker-Planck type, with 
quantum diffusion coefficient T> given by the square of the modulus of the recoil velocity 
matrix of the bound state system of D-particles and strings: 

9 

£> = caVa* X^l^l 2 (1-5) 

i=l 

where cq is a numerical constant and U l ab is the (renormalized) constant velocity matrix of 
a system of iV D-particles arising due to the D-particle recoil from the scattering of string 
states. This phenomenon is in fact characteristic of Liouville string theory, on which the 
above approach is based. Since the D-particle interactions distort their surrounding space- 
time, these non-linear structures may be thought of as describing short-distance quantum 
gravitational properties of the D-brane spacetime. Non-linear equations of motion for string 



field theories have been derived in other contexts in |16[]. From this nonlinear Schrodinger 
dynamics we shall describe a multitude of classes of solutions, using Galilean invariance of 
the D-brane dynamics which is a consequence of the corresponding logarithmic conformal al- 
gebra. We will show that bound state solutions do indeed exist for string couplings g s larger 
than a dynamically determined critical value. The effective bound state mass is likewise 
determined as an energetically induced shift of the static, BPS mass of the DO-branes. In 
fact, we shall find that there are essentially three different phases of the quantum dynamics 
in string coupling constant space. Below the critical string coupling the multiple D-brane 
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wavefunction is described by solitary waves, in agreement with the description of free D- 
branes as string theoretic solitons, while at the critical coupling the quantum dynamics 
is described by coherent Gaussian wavepackets which determine the appropriate quantum 
smearing of the multiple D-particle spacetime. These results are shown to be in agreement 



with the previous results concerning the structure of quantum spacetime JT^JTTJ . 

We close this section by summarizing some of the generic guidelines that we shall use in 
this paper for constructing a wavefunctional for the system of D-branes. We will use a field 
theoretic approach by identifying the Hartle-Hawking wavefunction 



;i.6) 



where Se is the effective Euclidean action. We shall discuss the extension to string theory 
and highlight the advantages and disadvantages of using this identification. We shall also 
identify the probability density with the genus expansion of an appropriate worldsheet o~- 
model: 

v = = [ Dx e ~ Sa[x] ( L7 ) 

genera 

The arguments in favour of this identification will be reality, and the occurence of statistical 
probability distribution factors which appear in the wormhole parameters after resummation 
of ( |1.7| ) over pinched genera. The Wilson-Polchinski worldsheet renormalization group flow, 
coming from the sum over genera as in ( |1.7[ ), yields a Fokker-Planck diffusion equation 

d t V = VV 2 V -V -J (1.8) 



where T> is the diffusion operator defined in (|1.5| ) in terms of (renormalized) recoil velocity 
matrices, and J is the associated probability current density. The equation ( |1.8|) will follow 
from the gradient flow property of the a-model /3-functions, which is also necessary for the 
Helmholtz conditons or equivalently for canonical quantization of the string moduli space. 

The knowledge of the Fokker-Planck equation ( |1.8| ) alone does not lead to an unambigu- 
ous construction of the wavefunction There are ambiguities associated with non-linear 
^-dependent phase transformations of the wavefunction: 

AT 7jA (^) = 7 log|*| + Aarg^ + #({>7 6 },t) (1.9) 

where t is the Liouville zero mode. Furthermore, \1/ is then necessarily determined by a 
non- linear wave equation if a diffusion coefficient T> is present, as will be the case in what 
follows. The non-linear Schrodinger equation has the form 
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if) X7 2 V 

ih d t y = H y + -v — v (1.10) 

where V — \E r t\E r is the probability density. This is a Galilean-invariant but time-reversal 
violating equation, exactly as expected from previous considerations of non-relativistic D- 
brane dynamics and Liouville string theory. Eq. (|1.10| ) will be the proposal in the following 
for the non-linear quantum dynamics of matrix D-branes (this was noted in passing in [|T7jj). 



II. QUANTUM MECHANICS ON MODULI SPACE 



In [|H]] it was shown how a description of non-abelian D-particle dynamics, based on 
canonical quantization of a a-model moduli space induced by the worldsheet genus expansion 
(i.e. the quantum string theory), yields quantum fluctuations of the string soliton collective 
coordinates and hence a microscopic derivation of spacetime uncertainty relations, as seen by 
short distance D-particle probes. In the following we will proceed to construct a wavefunction 
for the system of DO-branes which encodes the pertinent quantum dynamics. To start, in 
this section we shall clarify certain facts about wavefunctionals in non-critical string theories 



in general, completing the discussion put forward in 115]. 



A. Liouville-dressed Renormalization Group Flows 

Consider quite generally a non-critical string a-model, defined as a deformation of a 
conformal field theory with coupling constants {g 1 }- The worldsheet action is 

S a [x; {g 1 }] = S.[x] + J d 2 z g'V^x] (2.1) 

s 

where Vj are the deformation vertex operators and an implicit sum over repeated upper and 
lower indices is always understood. We assume that the deformation is relevant, so that the 
worldsheet theory must be dressed by two-dimensional quantum gravity in order to restore 
conformal invariance in the quantum string theory. The corresponding Liouville-dressed 
renormalized couplings {A 7 } satisfy the renormalization group equations 

X 1 + QX 1 = -f3\X) (2.2) 

where the dots denote differentiation with respect to the worldsheet zero mode of the Liou- 
ville field. Here Q is the square root of the running central charge deficit on moduli space 
and 

= tfX 1 + c r JK X J X K + . . . (2.3) 
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are the flat worldsheet /^-functions, expressed in terms of Liouville-dressed coupling con- 
stants. In (|2.3|), h 1 are the conformal dimensions and c 1 JK the operator product expansion 
coefficients of the vertex operators Vj. The minus sign in ( |2.2| ) owes to the fact that we con- 
fine our attention here to the case of central charge c > 25 (corresponding to supercritical 
bosonic or fermionic strings). 

Upon interpreting the Liouville zero mode as the target space time evolution parameter, 
eq. ( |2.2p is reminescent of the equation of motion for the inflaton field in inflationary 
cosmological models |18,n|. In the present case of course one has a collection of fields {g 1 }, 
but the analogy is nevertheless precise. The role of the Hubble constant H is played by the 
central charge deficit Q. The precise correspondence actually follows from the gradient flow 
property of the string a-model /3-functions for flat worldsheets: 

(3 1 = G IJ ^jC (2.4) 
dg J 

where C = Q 2 is the Zamolodchikov C-function which is associated with the generating 
functional for one-particle irreducible correlation functions [2D , and G IJ is the matrix inverse 
of the Zamolodchikov metric 

G IJ = 2\z\ 4 (V I (z,z)Vj(0,0)} (2.5) 



on the moduli space M({g 7 }) of cr-model couplings {g 1 }- Then the right-hand side of ( |2.2| ) 
also corresponds to the gradient of the potential V in inflationary models: 

+ 3*0=-^ (2.6) 

where is the inflaton field in a sufficiently homogeneous domain of the universe. 



B. The Hartle-Hawking Wavefunction 
In [|l~0| ) [r5|l it was shown, through the energy dependence of quantum uncertainties, that 



some sort of stochasticity characterizes non-critical Liouville string dynamics, implying that 
the analogy of eq. ( [2.2|) with the equations of motion in inflationary models should be made 



with those involving chaotic inflation ||19|| . Let us now briefly review the properties of these 
latter models. In such cases, the ground state wavefunction of the universe may be identified 

as 



21 



ip (a, <p) = exp -S E (a, (f>) (2.7) 

where Se is the Euclidean action for the scalar field a(r) and the inflaton scalar field 0(r) 
which satisfy the boundary conditions: 



S 



a(0) = a , 0(0) = 



(2.8) 



and r is the Euclidean time. 

To understand how eq. ( |2.7| ) comes about, we appeal to the Hartle-Hawking interpreta- 
tion JO]]. Consider the Green's function (x, t|0,t') of a particle which propagates from the 
spacetime point (0, t') to (x, t): 

(x, t\0, t'} = ]T il>l( x )M Q ) e iEn ^ = I Dx e^'W* (2.9) 

where {ip n } is the complete set of energy eigenstates with energy eigenvalues E n > (the 
sum in (|2.9|) should be replaced by an appropriate integration in the case of a continuous 
spectrum). To obtain an expression for the ground state wavefunction, we make a Wick 
rotation t = —ir, and take the limit r — > — oo to recover the initial state. Then in the 
summation over energy eigenvalues in (|2.9|) , only the ground state (n = 0) term survives 
if Eq = 0. The corresponding path integral representation becomes / Dx e~ SE ^ x \ and one 
obtains eq. (|2.7| ) in the semi-classical approximation. 

For inflationary models which are based on the de Sitter spaces dS^ with 

a(r) = ^(0) cos/t(0)r (2.10) 

one has 



and hence 



s ^> = -w(0) (2 - n; 



* (a ^ ) = exp T6FW (2 - 12) 



Thus the probability density for finding the universe in a state with = const., a = k 1 (<p) 



' is 



8irV((j>) 



The distribution function ( 2.131 ) has a sharp maximum as V(0) — > 0. For inflationary models 



this is a bad feature, because it diminishes the possibility of finding the universe in a state 
with a large field and thereby having a long stage for inflation. However, from the point 
of view of Liouville string theory, the result Q2.13 ), if indeed valid, implies that the critical 



string theory (since V oc Q 2 there) is a favourable situation statistically, and hence any 
consideration (such as those in jn|) made in the neighbourhood of a fixed point of the 
renormalization group flow on the moduli space of running coupling constants is justified. 
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C. Moduli Space Wavefunctionals 



Let us now proceed to discuss the possibility of finding a Schrodinger wave equation for 
the D-particle wavefunction. The identification (277) in the inflationary case needs some 



careful verification in the case of the topological expansion of the worldsheet a-model ( j2.1fl . 
In Liouville string theory, the genus expansion of the partition function may be identified 
|T5[ with the wavefunctional of non-critical string theory in the moduli space of coupling 
constants {g 1 }: 

#({/})= Yl /^e-W"=e-« (2.14) 

genera 

where 

oo 

W}] = EG^~ 2 -W}] (2-15) 

h=0 

is the effective target space action functional of the non-critical string theory. The sum 
on the right-hand side of ( |2.15 ) is over all worldsheet genera, which sums up the one- 



particle irreducible connected worldsheet amplitudes Th with h handles. The gradient flow 



property (|2.4j) of the /3-functions ensures [10,15] that the Helmholtz conditions for canonical 



quantization are satisfied, which is consistent with the existence of an off-shell action jF[{g 7 }]. 
In that case, the effective Lagrangian on moduli space whose equations of motion coincide 
with the renormalization group equations (|2.2|) is given by JTD 



C M {t) = -(3 I G IJ (3 J (2.16) 

and it coincides with the Zamolodchikov C-function. The semi-classical wavefunction deter- 
mined by ( |2.14j ) is thereby determined by the action C[A] regarded as an effective action on 



the space of two-dimensional renormalizable field theories. Thus the probability density is 
^[{fl^}] = e^ 2 ^ 9 '^, which implies that the minimization of .^[{g 7 }] yields a maximization 
of provided that the effective action is positive-definite. This is an ideal situation, 

since then the minimization of .Ff-fg 7 }], in the sense of solutions of the equations SJ-'/Sg 1 = 0, 
corresponds to the conformally-invariant fixed point of the cr-model moduli space, thereby 
justifying the analysis in a neighbourhood of a fixed point. 

However, the identification ( |2.14j ) is not the only possibility in non-critical string theory, 
as will be discussed below, in particular in connection with the Schrodinger dynamics of 
DO-branes. The main point is that upon taking the topological expansion in Liouville string 
theory, the couplings g 1 become quantized in such a way so that 

£ ' f Dx e-^te^ = / Da 1 e"^^ / Dx e'^W^}] {2 U) 

otiord J J J 



genera 



M({aO) 
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where the prime on the sum means that the genus expansion is truncated to a sum over 
pinched annuli of infinitesimal strip size, Sf'[x; {g 1 }} is the tree-level (disc or sphere) action 
for the u-model, and a 1 are worldsheet wormhole parameters on the moduli space M({g 7 }) 
of the two-dimensional quantum field theory. The Gaussian spread in the a 1 in ( |2.17| ) can 
be interpreted as a probability distribution characterizing the statistical fluctuations of the 
coupling constants g 1 . The width V is proportional to the logarithmic modular divergences 
on the pinched annuli, which may be identified with the short- distance infinities log A at 
tree-level [[Hj (A is the worldsheet ultraviolet cutoff scale). The result (|2.17| ) suggests that 
one may directly identify the genus expansion of the worldsheet partition function as the 
probability density 

*({/},*)| 2 =?W},t) (2-18) 

for finding non-critical strings in the moduli space configuration {g 1 } at Liouville time t (the 
worldsheet zero mode of the Liouville field). In this way one has a natural explanation for 



the reality of eq. ( [2. 17| ) on Euclidean worldsheets. If the identification of the genera summed 
partition function with the probability density holds, i.e. with the square of the wavefunction 
^{{g 1 }^) rather than the wavefunctional itself, then one may obtain a temporal evolution 
equation for ( |2.18| ) using the Wilson-Polchinski renormalization group equation on the string 
worldsheet WM. This will be described in section IV. 

One may argue formally in favour of the above identification in the case of Liouville 
strings, within a world-sheet formalism, by noting [2~2 that the conventional interpretation of 
the Liouville (world-sheet) correlators as target-space S-matrix elements breaks down upon 
the interpretation of the Liouville zero-mode as target time. Instead, the only well-defined 
concept in such a case is the non-factorizable ^-matrix, which acts on target-space density 
matrices rather than state vectors. This in turn implies that the corresponding world-sheet 
partition function, summed over topologies, which in the case of critical strings would be 
the generating functional of such S'-matrix elements in target space, should be identified 
with the probability density in the moduli space of the non-critical strings ( |2.18| ). In the 
Appendix we review this approach by focusing on those aspects of the formalism that 
are most relevant to our pruposes here. As we shall discuss there, the above identification 
follows from specific properties of the Liouville string formalism. 

Notice that if one interprets the topological expansion of the worldsheet partition function 
as the probability density for the non-critical string configuration {g 1 }, then the simple 
argument leading to eq. ( |2.7| ) is not valid here. In such a situation the action in eq. (f2.9|), 
which refers to the string moduli space, is not the same as the effective target space action 
jF[{(j( / }], but rather something different, corresponding to the phase of the wavefunctional 
^({g 1 }^) whose probability density (|2.18| ) corresponds to the genera summed worldsheet 
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partition function. This is not necessarily a bad feature, as we shall see, although in most 
treatments the target space effective action ^[{p 7 }] i s identified with the moduli space action 
upon identification of the Liouville zero mode (i.e. the local worldsheet renormalization 
group scale) with target time. For this, we observe that the statistical interpretation of the 
resummed worldsheet partition function is compatible with the interpretation in [K| of the 
Gaussian wormhole parameter distribution function in eq. ( J2.17 ) as being responsible for 
the quantum uncertainties of D-branes. This follows trivially from the fact that 



-2.F({ ff '},t) 



(2.19) 



Then, any correlation function may be written as 
(V h ---V In ) = J Dg 1 ^{{g 1 }^ V h ---V In 

M({gl}) 

J Dg 1 J Da 1 e^'^ J Dx e'^^ 1 ^ V h ■ ■ ■ V In (2.20) 



gl}) M({gl}) 

which using eq. ( |2.19| ) gives the connection between the two probability distributions. 



III. MATRIX D-BRANE DYNAMICS 



In this section we shall briefly review the worldsheet description of [T(| for matrix D0- 
brane dynamics. The partition function is given by || 



Z[A ,Y}= J D(jl(x,£,£) exp 



J d 2 z v, &c" dx v + ^j dr x t {r) d a x\r) 
s as 



(3.1) 



where 



W[x, £, £] = exp ig s f dr ^^(r) d T x\r) + — ; U^Yf (x u )&(r) d a x\r) J (3.2) 

is the deformation action of the free cx-model in ( |3.1|) . Here the indices fi = 0, 1, . . . , 9 and 
i — 1, ... ,9 label spacetime and spatial directions of the target space, which we assume has 
a flat metric rj^. The functional integration measure in ( |3.1| ) is given by 



D^{x,li) = Dx^ D£D£ exp 



J2 ( fdr Ur) d T Ur) + ^(0)^(0) 



N 



Ee-(o)e.(i) 

(3.3) 



a=l 
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The complex auxilliary fields ^ a (r) and £ q (t), a = 1, . . . , N, transform in the fundamental 
representation of the brane gauge group, and they live on the boundary of the worldsheet 
X which at tree-level is a disc whose boundary is a circle <9£ with periodic longitudinal 
coordinate r G [0, 1] and normal coordinate a G R. They have the propagator (Ca( T )6&( r/ )) — 
S a b @(t' — t), where denotes the usual step function. The integration over the auxilliary 
fields with the measure ( |3.3| ) therefore turns Q3.2| ) into a path-ordered exponential functional 
of the fields x which is the T-dual of the usual Wilson loop operator for the ten-dimensional 
gauge field (A , —-^t^Y 1 ) dimensionally reduced to the D-particle worldlines. In this picture, 
A is thought of as a gauge field living on the brane worldline, while Yf^, a = 1, . . ., N, 
are the transverse coordinates of the iV D-particles and Y°" b , a ^ b, of the short open 
string excitations connecting them. We shall subtract out the center of mass motion of the 
assembly of iV D-branes and assume that Yi G su(N). We shall also use S , £/(iV)-invariance 
of the theory ( [3.1|) to select the temporal gauge A = 0. 

The action in ( |3.1| ) may be formally identified with the deformed conformal field theory 
( j2.ip by taking the couplings g 1 ~ and introducing the one-parameter family of bare 
matrix-valued vertex operators 

V: b (x; r) = ^ d a x\r) &(t)&(t) (3.4) 

This means that there is a one-parameter family of Dirichlet boundary conditions for the 
fundamental string fields x l on <9£, labelled by r G [0, 1] and the configuration fields 

y l {x -r)=Ur)Yt\x\r))Ur) (3.5) 

Instead of being forced to sit on a unique hypersurface as in the case of a single D-brane, in 
the non-abelian case there is an infinite set of hypersurfaces on which the string endpoints 
are situated. In this sense the coordinates ( |3.5|) may be thought of as an "abelianization" 
of the non-abelian D-particle coordinate fields Y^ h . 

To describe the non-relativistic dynamics of heavy D-particles, the natural choice is 
to take the couplings to correspond to the Galilean boosted configurations Y°- b (x Q ) = 

yah + Jjab x 0^ w ^ eIf} JJ. t h e 

non-relativisitic velocity matrix. However, logarithmic modular 
divergences appear in matter field amplitudes at higher genera when the string propagator 
Lo is computed with Dirichlet boundary conditions. These modular divergences are can- 
celled by adding logarithmic recoil operators [0,0 t° the matrix a-model action in ( |3.1|) . 



From a physical point of view, if one is to use low-energy probes to observe short- distance 
spacetime structure, such as a generalized Heisenberg microscope, then one needs to con- 
sider the scattering of string matter off the assembly of D-particles. For the Galilean-boosted 
multiple D-particle system, the recoil is described by taking the deformation of the cx-model 
action in (|3.1| ) to be of the form 
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Y t ab (x°) = VtfY^C^x ) + UfD e (x°) = (VdeYf b + Ufx°) e (x°) (3.6) 

where 

C e (x°) = e& e (x°) , D e (x°) = x° & e (x°) (3.7) 

and 

1 + °° A 

e e ( x °) = — ( e^° (3.8) 

K J 2m J q-ie K J 

— oo 

is the regulated step function whose e — > + limit is the usual step function. The operators 
( |3.7| ) have non-vanishing matrix elements between different string states and therefore de- 
scribe the appropriate change of quantum state of the D-brane background. They can be 
thought of as describing the recoil of the assembly of D-particles in an impulse approxima- 
tion, in which it starts moving as a whole only at time x° = 0. The collection of constant 
matrices {Y* b , U 3 cd } now form the set of coupling constants {g 1 } for the worldsheet a-model 

The recoil operators (|3.7|) possess a very important property. They lead to a deformation 



of the free cx-model action in (|3.1|) which is not conformally-invariant, but rather defines a 
logarithmic conformal field theory [p4| . Such a quantum field theory contains logarithmic 
scaling violations in its correlation functions on the worldsheet, which can be seen in the 



present case by computing the pair correlators of the fields (3/7) [23 

[C e (z)C e (0)) = 

;C e (z)£> e (0) 



b 



D £ (z)D e (0)) = ^logz (3.9) 



where 



lei 2 a' 



K = (3.10) 

is the conformal dimension of the recoil operators. The constant b is fixed by the leading 
logarithmic divergence of the conformal blocks of the theory. Note that ( |3.10| ) vanishes as 
e — > 0, so that the logarithmic worldsheet divergences in ( |3.9| ) cancel the modular annulus di- 
vergences mentioned above. An essential ingredient for this cancellation is the identification 



\ = -2a' log A (3.11) 



14 



which relates the target space regularization parameter e to the worldsheet ultraviolet cutoff 
scale A. 

Logarithmic conformal field theories are characterized by the fact that their Virasoro 
generator Lq is not diagonalizable, but rather admits a Jordan cell structure. Here the 
operators ( |3.7| ) form the basis of a 2 x 2 Jordan block and they appear in the spectrum 
of the two-dimensional quantum field theory as a consequence of the zero modes that arise 
from the breaking of the target space translation symmetry by the topological defects. The 
mixing between C and D under a conformal transformation of the worldsheet can be seen 
explicitly by considering a scale transformation 



A -> A' = Ae" 



(3.12) 



Using ( |3.11| ) it follows that the operators ( |3.7| ) are changed according to D' e = D e + tya'C t 
and C' e = C e . Thus in order to maintain scale-invariance of the theory (|3.1| ) the coupling 
constants must transform under ( |3.12|) as [23 , ^5 Y n = Y % + U l t and U n = U l , which are 
just the Galilean transformation laws for the positions Y % and velocities U i . Thus a scale 
transformation of the worldsheet is equivalent to a Galilean transformation of the moduli 
space of cr-model couplings, with the parameter e~ 2 identified with the time evolution pa- 
rameter t = — v / a'logA. The corresponding /5-functions for the worldsheet renormalization 
group flow are 



A» S f = M5 + 



dt 



KUi 



(3.13) 



and they generate the Galilean group G(9) n2 in nine-dimensions. 
The associated Zamolodchikov metric 



G% cd = 2NA 2 (V: b (x; 0) V> d (x; 0) 



(3.14) 

can be evaluated to leading orders in cr-model perturbation theory using the logarithmic 

and the propagator of the auxilliary fields to give [ID 



conformal algebra 



4o 2 

pi] _ 

" ab;cd i 

a' 



In ® In + 36 { In ® + f/if/J ) 

+ 1/* ® u j + u j ®u' + (fru* + WU') ®I N \] +0 (f s ) 

V ' J db:ca V ' 



(3.15) 



where ijv is the identity operator of SU (N) and we have introduced the renormalized cou- 
pling constants 



g s = g s / Va'e 



U l = U l /\fa'e 



(3.16) 
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From the renormalization group equations ( |3.13|) it follows that the renormalized velocity 
operator in target space is truly marginal, 







(3.17) 



which ensures uniform motion of the D-branes. It can also be shown that the renormalized 
string coupling g s is time-independent |T(J . If we further define the position renormalization 



(3.18) 



then the /3-function equations ( ^.13| ) coincide with the Galilean equations of motion of the 
D-particles, i.e. 



dY i 
~~dT 



IP 



(3.19) 



Note that the Zamolodchikov metric ( |3.15| ) is a complicated function of the D-brane dynam- 
ical parameters, and as such it represents the appropriate effective target space geometry of 
the D-particles. The moduli space Lagrangian (|2.16|) is then readily seen to coincide with 
the expansion to 0{g^) of the symmetrized form of the non-abelian Born-Infeld action for 
the D-brane dynamics 0, 



^nbi = 7= — - tr Sym /det [rj^ I N + 2na'g 2 s F^] 
V Ina' g s V ^ v 

where tr denotes the trace in the fundamental representation of SU(N), 

1 



(3.20) 



Sym(M 1 ,...,M„ 



in 



(3.21) 



is the symmetrized matrix product and the components of the dimensionally reduced field 
strength tensor are given by 



1 dY { 



2rca' dt 



Y Y 

* ' 3 



(3.22) 



IV. EVOLUTION EQUATION FOR THE PROBABILITY DISTRIBUTION 

In this section we will derive the temporal evolution equation for the probability density 
V({g J },t) following the identification of time with a worldsheet renormalization group scale 
(i.e. the Liouville zero mode). The basic identity is the Wilson-Polchinski equation for the 
case of the worldsheet action (|2.1|) which reads |14[ 
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dZ 
d log A 

Dx» e - s ° lx ' i9l}] 



d z id w 



dS: 



hit 



Slog A 
d 



d log A 



G(z-w) 



5 2 S; 



hit 



+ 



hit 



5x^(z)5x ll (w) 5x^(z) Sx^w) 



(4.1) 



and it is the requirement of conformal invariance of the quantum string theory. Here S- m t 
S a — Z is the partition function of the cx-model, and 



G(z-w) = {° xP(z)xr(w)l) m 



(4.2) 



is the two-point function computed with respect to the conformal field theory action S* [x] . 
The basic assumption in arriving at eq. ( [4.1| ) is that the ultra-violet cutoff A on the string 
worldsheet appears explicitly only in the propagator G(z — w), as can always be arranged 
by an appropriate regularization |L4] . 



Henceforth we shall concentrate on the specific case of interest of a system of N interact- 
ing D-particles. Then, upon summing up over pinched genera, there are extra logarithmic 



divergences in the Green's function ( [4.2[) coming from pinched annulus diagrams, which may 
be removed by the introduction of logarithmic recoil operators, as explained in the previous 
section. Using primes to denote the result of resumming the topological expansion over 
pinched genera, we then have that 







Slog A 



G(z-w)' 







d log A 



d 



genera 



d log A 



where the correlator {-)- m t includes the disc and recoil interaction contributions. Subtracting 
the disc A-dependence in normal ordering, the remaining dependence on the worldsheet 
cutoff comes from the two-point functions of the logarithmic recoil operators, giving terms 
of the form 







Slog A 



° o x»(z)x,(w) ° (a C cC e (z)C e (w) + a OD C e (z)D e (w) + a DD D e (z)D e (w)))^ (4.4) 



The leading divergence comes from the correlation function (D t (z)D e (w))* ~ logA, which 
follows upon the identification (|3.11|) . Thus we may write 







9 N 



G(z-w)~c G (a'f\og\z-w\Y^ E 1^ 

UI °S 1V i=l a ,b=l 



i 1 2 
ab I 



(4.5) 



where cq > is a numerical coefficient whose precise value is not important, and we have 
used the fact that U l G su{N). 
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Next, we observe that in the case of D-particles the second term in eq. ( |4.1| ) becomes 

9 N 

J D»(x,l£) e~ s « jdr jdr' (a') 2 c G J2 £ log[2 - 2 cos(t - r')] 



x 



as dT, i=1 a < b=1 



Sx^^Sx^t') &c^(t) &c M (r' 
where the interaction Lagrangian is given by 



Sm = ^f dT ^{t) Ur) Yf(x°) &(t) (4.7) 

In the case of a system of recoiling DO-branes, the cx-model couplings in eq. ( (4.7|j are given 
by ( |3.6| ) with the abelianized couplings fl3.5| ) of Y^ 6 viewed as the boundary values for the 
open string embedding fields x % (r) on the D-brane. This means that the fields x*(r) are 
simply identified with ^a^Y^^T). All the non-trivial dependence comes from the x° field 
which obeys Neumann boundary conditions and is not constant on the boundary of E. Then 
we may write 



(4.6) 



8x^(t)5x^t') 5x»{t) Sx^t') 

= V^ int + (V^ int ) 2 + (^) 2 UfUfUr)Ur)Ur')Ur) 

x a (7 x i (T)a (J x• ? '(r , )e £ (s (r))e £ (x (r , )) (4.8) 

where z/i denotes the constant abelianized zero modes of x l (r) on <9£. Here we have used 
the fact that terms of the form x°5(x°) and Q(x°)5(x°) vanish with the regularization ( |3.8|) 
Toll . The terms involving d cr x l (T)d cr x^ (r') will average out to yield terms of the form 



\K b \ 2 (e e (x°(r))e e (x (r')))^ = a'\U* ab \ 2 (C e (r)C e (/)>, ~ 0(e 2 ) (4.9) 

where we have used the logarithmic conformal algebra. At leading orders, these terms vanish, 
but we shall see the importance of such sub-leading terms later on. 
Using the Dirichlet correlator 

^(r)fe(rO) =- 36?rV „ (4.10) 
' * 1 — cos (r — r ) 



we find that the boundary integrations in eq. (|4.6|) are of the form [10 

, r , , log[2-2cos(r-r / )] . . 

dr i dr' — — ^ ~ log A (4.11) 

J 1 — cos(r — r ) 

<9£ dT, y ' 

which has the effect of renormalizing the velocity matrix U l ab — > U l ab . Thus, ignoring the 
0(e 2 ) terms for the moment, we find that the remaining terms in the Wilson- Polchinski 
renormalization group equation (|4.1| ) yield a diffusion term for the probability density: 
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9 N 

d t V[Y,U;t} = c G V^Yl E K b \ 2 VlV[Y,U;t] + 0(e 2 ) (4.12) 

j=l a,6=l 

This equation is of the Fokker-Planck type, with diffusion coefficient 

9 N 

v= Cg v^y: e \K b \ 2 (4.i3) 

i=l a, 6=1 

coming from the quantum recoil of the assembly of D-particles. The diffusion disappears 
when there is no recoil. Note that ( 4.13|) naturally incorporates the short-distance quantum 
gravitational smearings for the open string interactions (compare with eq. (|1.2fj ), and it 
arises as an abelianized velocity for the constant auxilliary field configuration $, a (r) = £ (t) = 
1 Va = 1,...,N. 

The evolution equation ( |4.12j ) should be thought of as a modification of the usual conti- 
nuity equation for the probability density. Indeed, as we will now show, the 0(e 2 ) terms in 
eq. ( 4.12j ) coming from (|4.9|) are of the form — V w j7i, where 

is the probability current density. Here 

m = -j=r , = 4<? s (4.15) 

Va'g s 

are, respectively, the BPS mass of the D-particles and the moduli space "Planck constant" . [] 
For this, we note first of all that such terms should generically come in the form 

- V w Ji = (Vl arg*) V-^ (V w arg*) V m V (4.16) 



The identification (4.15) of Planck's constant in the D-particle quantum mechanics on moduli 
space with the string coupling constant is actually not unique in the present context of considering 
only the exchange of strings between D-particles. As discussed in |fLO| , the most general relation, 
compatible with the logarithmic conformal algebra, involves an arbitrary exponent x through 
= 4(g s ) 1+x ^ 2 . The exponent x arises from specific mechanisms for the cancellation of modular 
divergences on pinched annular surfaces by appropriate world-sheet short-distance infinities at 
lower genera. The only restriction imposed on x ls that it be positive definite. As shown in [1C], 
the standard kinematical properties of D-particles are reproduced by the choice x = | . A choice of 
X 7^ seems more natural from the point of view that modular divergences should be suppressed 
for weakly interacting strings. However, in the present case, we assume for simplicity the value 
X = 0, which yields the standard string smearing \fa! for the minimum length uncertainty. The 
incorporation of an arbitrary x ^ in the formalism is straightforward and would not affect the 
qualitative properties of the following results. 
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The second term in eq. ( |4.16|) , upon identification of the probability density V with the 
genera resummed partition function on the string worldsheet, is proportional to the world- 
sheet renormalization group /3-function, given the gradient flow property (|2.4j) of the string 



effective action [BOfl, so that 



(4.17) 



which is to be understood in terms of abelianized quantities. In the present case the renor- 
malization group equations are given by (|3.17|) and (|3.19|) and, since the couplings Uf' 
are truly marginal, we are left in (|4.17|) with only a Zamolodchikov metric contribution 
Gcc = 2NA A {C t (T)C t (r)) 3f (Note that here one should use suitably normalized correlators 
(•)* which yield the behaviour (|4.17p ). From the logarithmic conformal algebra it therefore 
follows that a term with the structure of the second piece in eq. (|4.16 ) is hidden in the 
contributions ( |4.9|) which were dropped as being subleading in e. Furthermore, from ( |3.15| ), 
( CT ), ( EL5D , (gjSD and O it follows that to leading orders 



V yi arg ^ 




(4.18) 



where U; 



dt/i/dt is the worldsheet zero mode of the abelianized, renormalized velocity 

2 



operator. It then follows that to leading orders we have V 2 . arg \P = 0.0 

Thus, keeping the subleading terms in the target space regularization parameter e leads 
to the complete Fokker-Planck equation for the probability density V = ^>^: 



d t V[Y, U- 1] = -V Vi Ji\V, U;t}+V VlV[Y, U; t] 



(4.19) 



where Ji is the probability current density (|4.14|) , with \1/ the wavefunctional for the system 
of D-branes: 



*[Y,[/;t] = nexp 



i=l 



1C G Vi 



E tr l^l 



W,U;t] 



(4.20) 



2 Noncommutative position dependent terms arising from commutators \Yi, Yj] appear only at two- 
loop order in cr-model perturbation theory 0]. An interesting extension of the present analysis 
would be to generalize the results to include these higher-order terms into the quantum dynamics. 
However, given that the pertinent equations involve only the abelianized coordinates ([Tl|), we do 
not expect the inclusion of such terms to affect the ensuing qualitative conclusions. The effect 
of the noncommutativity is to render the quantum wave equation for the system of D-particles 
non-linear, through the recoil-induced diffusion from the multi-brane interactions, as we discuss in 
the subsequent sections (for a single brane one would obtain a free wave equation governing the 
quantum dynamics). 
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Such quantum diffusion is characteristic of all Liouville string theories [|T^,|T7],^(J . The re- 
sulting quantum dynamics, including the quantum diffusion which arises from the D-brane 
recoil, is described by the Schrodinger wave equation which corresponds to this Fokker- 
Planck equation. This equation is analysed in detail in the next section. 



V. NON-LINEAR SCHRODINGER WAVE EQUATIONS 



Given the Fokker-Planck equation ( |4.19|) , there is no unique solution for the wavefunction 



as we discuss below, and the resulting Schrodinger wave equation is necessarily non-linear, 



due to the diffusion term [27,28]. Consider the quantum mechanical system with diffusion 



which is described by the Fokker-Planck equation ( }4.19|) for the probability density V = Wty. 



In [27] it was shown that, by imposing diffeomorphism invariance in the space y E M and 
representing the symmetry through the infinite-dimensional kinematical symmetry algebra 
C°°(M) -B Vect(M), one may arrive at the following non-linear Schrodinger wave equation: 

in M -^- = Wo* + i/(*)* (5-1) 
where Ti,Q is the linear Hamiltonian operator 

H = -^Vl + V u (y,u;t) (5.2) 

and 

Here Vm(v, u; t) is the interaction potential on moduli space and the real continuous quantum 
number D in (|4.13| ) is the classification parameter of the unitarily inequivalent diffeomor- 



phism group representations. Other models which have more than one type of diffusion 
coefficient can be found in [[27],^ . 



A crucial point [p£| is that there exist non-linear phase transformations of the wave- 
function \& (known as quantum mechanical "gauge transformations") which leave invariant 
appropriate families of non-linear Schrodinger equations, and also the probability density V. 
Such transformations do not affect any physical observables of the system. This implies that 
the choice of \I/ is ambiguous, once a density V is found as a solution of eq. ( |4.19| ) on the 
collective coordinate space {Yj afe } of the D-branes. An important ingredient in finding such 
transformations is the assumption p8| , |29[] that all measurements of quantum mechanical sys- 
tems can be made so as to reduce eventually to position and time measurements. Because of 
this possibility, a theory formulated in terms of position measurements is complete enough 
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in principle to describe all quantum phenomena. This point of view is certainly met by 
the D-brane moduli space, whereby the wavefunctional depends only on the couplings {g 1 } 
and not on the conjugate momenta pi = —ihy^d/dg 1 . The group of non-linear gauge trans- 
formations acts on each leaf in a foliation of a family of non-linear Schrodinger equations, 
such that the two-dimensional leaves of the foliation consist of sets of equivalent quantum 
mechanical evolution equations. 

It follows that then one can perform the following local, two-parameter projective gauge 
transformation of the wavefunction p8fl: 



= N 7 ,a(*) = 1*1 exp(ry log |^| + iAarg^) (5.4) 

under which the probability density is invariant, but the probability current transforms as 

Jl = \Ji + \v yi V (5.5) 

Here 7(t) and X(t) ^ are some real- valued time-dependent functions. The collection of all 
non-linear transformations iV 7i A obeys the multiplication law of the one-dimensional affine 
Lie group Aff(l). Under (5T4") there are families of non-linear Schrodinger equations that 



are closed (in the sense of "gauge closure"). A generic form of such a family, to which the 
non-linear Schrodinger equation ( |5.1| ) belongs, is 

*5r = ^-Wo* + w2-R3[*]* + /ii-Ri[*]*+ (a* 2 - R2M * 
at Tim 

+ (/i 3 + vi)R 3 [V] * + ^R^>] * + Us + \ vi) Rs[*] * 



i=l,2 3=1 ft M 

where z/j,/ij are real-valued coefficients which are related to diffusion coefficients T> and V 
by 

2m 



u 2 = -V 
2 

/i! = c{D' 

A*2 = ~i h c 2 £> 

4m 

/^3 = 77" + 

2m 

A*6 = ^ + c 5 ^' (5.7) 
8m 

and -Rj[*] are non-linear homogeneous functionals of degree which are defined by 
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R = m VyiJi 



R 2 
R, 

i?4 



v 

tl M r 



-p2 



In eq. ( |5.7|) the Cj are constants, while in eq. (|5.8|) the probability current density is given 
by ( pl|) with V = 

The gauge group Aff(l) acts on the parameter space of the family ( |5.6|) . Some members 
of this family are thereby linearizable to an ordinary Schrodinger wave equation under the 
action of (|5.4|) . These are the members for which there exists a specific relation between T> 
and T>' ||28| , and for which Ehrenfest's theorem of quantum mechanics receives no dissipative 
corrections. The quantum mechanics of D-particles is not of this type, given that there 
is definite diffusion, dissipation and thus time irreversability. However, as discussed in 
T0| , |23|p5| , one needs to also maintain Galilean invariance, which is a property originating 



from the logarithmic conformal algebra of the recoil operators. As described in [p^j , there is a 
class of non-linear Schrodinger wave equations which is Galilean invariant but which violates 
time-reversal symmetry. For this, it is useful to first construct a parameter set of equations 
of the form ( |5.7D which remain invariant under the gauge transformations (|5.4j ). We may 
describe the parameter family of equations ( |5.6|) in terms of orbits of Aff(l) by regarding 
7 = 2m/ii and A = 2mv\ as the group parameters of an Aff(l) gauge transformation Q5.4D . 
Then the remaining five parameters in ( |5.7|) are taken to be the functionally-independent 
parameters rjj, j = 1, . . . , 5, which are invariant under Aff(l) and are defined by 

1 

m = »2 - - mi 

^3 

V3 = — 
V\ 

^3 

774 = - Ml — 
V\ 

V5 = v\Vh ~ V2^a + fa) 2 — (5.9) 

^1 



A detailed discussion of the corresponding physical observables is given in [28|. For our 



purposes, we simply select the following relevant property of the non-linear Schrodinger 
equation based on the parameter set (p79|). 
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Consider the effect of time-reversal on the non-linear Schrodinger wave equation. Setting 
t — > — t is equivalent to introducing the following new set of coefficients: 

{h>i) T = -Vi £=1,2 
(Vj) T = ~H j = 1, . . . , 5 

(Vuf = -V M (5.10) 

where the superscript T denotes the time-reversal transformation. It is straightforward 
to show that, in terms of the rj/s, there is time-reversal invariance in the non-linear 



Schrodinger equation if the two parameters rji and 774 are both non-vanishing. On the other 



hand, a straightforward calculation also shows p8| that Galilean invariance sets 774 = 0, 
thereby implying that a family of non-linear Schrodinger wave equations which is invariant 
under G(9) but not time-reversal invariant indeed exists. For a single diffusion coefficient 



T> 7^ 0, as in the case ( |4.13|) of recoiling D-branes, one may set T>'cj = (corresponding to 



the Aff(l) gauge choice [i\ = 0) and thereby obtain the set of gauge invariant parameters: 

m = ~V 
' 2 

r)2 = 2a' g 4 s 

V3 = "I 
^4 = 

Vs = -a'g 4 s -\v 2 (5.11) 



The parameter set (|5.11|) breaks time-reversal invariance, as expected from the non-trivial 



entropy production and decoherence characterizing the worldsheet renormalization group ap- 
proach to target space time involving Liouville string theory |[T0|JT5| , |30|] . But it does preserve 



Galilean invariance, as is required by conformal invariance of the non-relativistic, recoiling 
system of D-particles. 

One may therefore propose that the Fokker-Planck equation for the probability density 
V on the moduli space of collective coordinates of a system of interacting D-branes implies 
a Schrodinger wave equation for the pertinent wavefunctional which is non-linear, Galilean- 
invariant and has a time arrow, corresponding to entropy production, and hence explicitly 
broken time-reversal invariance. The existence of a dissipation T> oc tr 1 0"* | 2 , due to the 
quantum recoil of the D-branes, implies that the Ehrenfest relations acquire extra dissipative 
terms for this family of non-linear Schrodinger equations. For example, one can immediately 
obtain the relations |27 
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where jr/ X,_0 ' ) is the undissipative current density Q4.14Q . Note that the fundamental renor- 
malization group equations Q3.13|) receive no corrections due to the dissipation. The existence 
of extra dissipation terms in Q5.12j) in the Ehrenfest relation for the momentum operator 
pi = —ihfyi V^, which are proportional to tr |?7*| 2 , may now be compared to the generalized 
Heisenberg uncertainty relations that were derived in | I0|| . These extra terms are determined 
by the total kinetic energy of the D-branes and their open string excitations, and they show 
how the recoil of the D-brane background produces quantum fluctuations of the classical 
spacetime dynamics. The relationship with quantum uncertainty relations will be discussed 
in the next section. 



VI. SOLUTIONS OF THE MATRIX D-BRANE WAVE EQUATION 



In this section we shall discuss the situation governing some of the solutions |27]|ri| of the 
class of non-linear Schrodinger equations ( |5.(j| ) which will be of interest to us in the context 
of the quantum states of a system of many D-particles. We may parametrize the generic 
wavefunctional ^ by an amplitude 6\(y,t) and a phase 62(11, t) as 



* = exp (6>i + id 2 ) 



(6.1) 



which, on account of (|5.6Q , satisfy the following coupled system of non-linear partial differ- 



ential equations [31 



= OA - 2i/ a V;.0! - i/iV^a - Au 2 (Vyfitf - 2v 1 V yi 6 1 V n 6 2 
= d t 6 2 + 2/x 2 V^i + /iiV^ 2 + 4(/i 2 + /i 5 ) (V y A) 2 

+ 2(/i! + ndVyftxVviO* + /i 3 (V Vi 6 2 f + t~~ V M (y, u; t) 



(6.2) 



This system can be equivalently expressed in terms of the gauge invariant set of parameters 
T]j defined in ( p.9| ) as 



dA = 2^V^i + — V 2 J 2 + A Vl (VyA) 



m 



^y i 6l^y i 6 2 



d t 6 2 = -krm) 2 V y .6\ - 8m (rj 2 + V5 + rjtfji - (Vi) Vs) (^yfii) 
- 2r j4 V Vi 1 V yi 2 - ^ (V yi 6 2 f + V M (y, u; t) 



(6.3) 
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where we have selected the gauge \i\ = 0. Notice that there are ambiguities in the phase 
function 9 2 in ( |6.1| ), which imply that the non-linear Schrodinger equation fl5.6| ) and the 
amplitude and phase equations ( |6.2[ ) are not completely equivalent. However, any solution 
of the phase and amplitude equation yields a solution of the non-linear Schrodinger equation 



fl5.6p . This is the commonly accepted point of view pi |, and the one which we adopt in this 
paper. 



A. Mass Superselection Rules 

One of the key properties of the matrix D-brane system that we have exploited extensively 
thus far is its Galilean invariance (or equivalently that the matrix a-model deformation 
defines a logarithmic conformal field theory). It is worthwhile to first mention some basic 
facts concerning Galilean invariant quantum field theories in the non-dissipative case {T> = 0) 



32fl . Consider two non-relativistic spinless fields 0i and (j> 2 of masses mi and m 2 with a 
quartic (pi(p 2 interaction mediated by a two-body potential V(x — y). The corresponding 
Haag expansions of their asymptotic fields are determined by Haag amplitudes fi(x — x a ,t — 
t a ; x — xj,t — tf), where a = 1, 2 labels the two fields and / labels the different bound states 
of the system. Then the function Fj(xi — x 2 ) defined through the equal time relation 

f - - ~\ rnixi + m 2 x 2 \ ,^ _> . . . 

//(xi - x 2 ; xi - £/) = 6 ■ Fi(xi - x 2 ) (6.4) 

V mi + m 2 l 



satisfies the usual Schrodinger wave equation 



n (— + — ) V ^-x' + V(xt - x 2 ) 
2 Vmi m 2 J 1 2 



Fj = -e/Fj (6.5) 



2 V?7ii m 2 , 

for the stationary bound state problem. This shows that the modified Haag amplitude is 
exactly the Schrodinger wavefunction of the bound state, with the reduced mass of the two 
particle system. 

The important property used in arriving at the result (|6.5| ) is the Galilean invariance 
of the system. The unitary projective representations of the Galilean group that arise in 
non-relativistic quantum mechanics cannot be reduced to vector representations, in contrast 
to the situation for most other physically relevant groups, such as the Poincare and Lorentz 
groups. Explicit mass parameters appear in the phase factors of these representations which 
lead to the Bargmann superselection rule |33|] that the sum of the masses appearing in the 



kinetic terms of the Hamiltonian must be conserved in every physical process. This implies, 
in particular, that the mass of any bound state must simply be the sum of the masses of 
its constituents. Note that since we are dealing with a non-relativistic theory, the energy 
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is not equal to the mass and indeed the energy of a two-particle bound state differs from 
the sum mi + m 2 by the binding energy. Since in the case of a composite system of N 
heavy D-particles, the non-linear Schrodinger equation also stems from the very specific 
Galilean invariance of the matrix D-brane dynamics, one might expect this superselection 
rule to carry over to the present diffusive situation. This would then imply that the mass 
appearing in the effective kinetic term should be Nm = Ng~ l /\fa' 1 and the number of 
colours N can be absorbed into the normalization of the mass, i.e. the iV-body DO-brane 
system behaves just like a single D-brane. However, this is not what we shall find below 
when we derive the extended Haag decomposition which takes into account the diffusion 
terms, i.e. of the quantum recoil of the D-particles. The massive Dirichlet string exchange 
which generates contact interactions among D-branes and gives rise to a non-linear quantum 
mechanical evolution equation must correctly be incorporated in order to distinguish a single 
brane from the non-abelian matrix D-branes. In the following we shall derive the modified 
superselection criteria that arise in the present context. 



B. Stationary Solutions 

The first important class of solutions of interest to us here are the stationary solutions 
which are defined as usual p7| by d{P = 0. This is equivalent to assuming d t Q\ = 0. 
Consider the free case whereby the moduli space interaction potential is ignored, Vm = 0. 
Then, setting 9i(y,t) = const., one obtains from the amplitude equation in ( |6.3|) the plane 
wave solutions ^(y,t) = fy e\"^~ u ^f ^ w jth a modified dispersion relation obtained from 
the resulting phase equation in (|6.3|) : 



d t e 2 {y,t) =u(k) = l^ + cs&jk 2 (6.6) 

Notice that for the Galilean invariant but time- reversal violating case (|5.11|) , with T>' = 0, 
the dispersion relation ( |6.6|) is that of a non-relativistic massive free particle. 

The fact that such plane wave solutions exist is an important feature of the present 
approach to the description of the effective dynamics of a system of many D-particles. As 
discussed in sections IV and V, the pertinent wave equation of a system of N D-particles, 
interacting via the exchange of strings, is written down in terms of the abelianized positions 
yi in ( |3.5| ), rather than an iV-body quantum mechanics that one would expect from the 
effective Yang-Mills reduction arising in the standard low-energy description of multiple D- 
particle dynamics This seems to tie in nicely with the heuristic "fat brane" picture of the 



multiple D-brane system advocated in 153], according to which one treats the assembly of 



D-particles as a single one with string interactions now encoded in the diffusion term ( |4.13|) . 
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Had one instead started with the Yang-Mills matrix quantum mechanics from the onset, 
then one would have obtained a linear Schrodinger equation with interactions taken care of 
by the Yang-Mills potential. This is what had been studied in the literature so far ||, with 
no definite, general conclusions about the bound state problem for strings and D-particles 
(see [|J for recent developments). In the present approach, we effectively reduce the quantum 
mechanics to a single body problem (the dynamics of the "fat brane" P3j), at the cost of 
obtaining a non-linear wave equation. This explains quite naturally the dispersion relation 
( p.6|) in the case where string interactions between the D-branes are ignored. Note also that 
the energy is proportional to g 2 s and thereby depends on the open string excitations between 
the constituent D-particles. 

A non-trivial issue concerns the nature of the effective moduli space potential Vm(v, u; t) 
among the constituent D-particles of the "fat D-brane". In the context of the cr-model 
approach to the system of recoiling D-particles |IJ , the effective potential VJvj among the D- 
particles is determined by the Zamolodchikov C-function for the deformed cr-model at hand, 
which according to the discussion in and in section III, is given by the Born-Infeld La- 
grangian ( |3.20| ). In the context of the (abelianized) "fat brane" picture, therefore, one would 
expect that the relevant Born-Infeld effective action can be reduced to the one containing 



only velocity terms of the form tryl — \U % \ 2 , given that the coordinate-dependent commu- 
tator terms in ( |3.22[ ) may be neglected to a first-order approximation, i.e. to leading orders 
in cr-model perturbation theory |]10[ . Such a coordinate independent potential may then be 
shifted to zero, and hence the free particle case Vm = discussed above corresponds to the 
Galilean invariant multiple D-particle system in a first-order approximation. This is similar 
in spirit (but not identical) to a discussion of the ground-state properties of the multiple 
D-particle system, given that the absence of the non-commutative commutator terms in the 
Born-Infeld action (|3.20|) also guarantees unbroken supersymmetry in the standard effective 
Yang-Mills target space description of multiple D-particle systems |35|] . This commutative 
limit is a saddle-point solution of the full U (N) quantum gauge theory. 

However, the free case is not a satisfactory approximation if one wishes to include quan- 
tum fluctuations about such ground states, i.e. to properly incorporate string interactions 
between the D-particles which may enable them to form a bound state. Such a case would 
correspond to off-shell non-conformal backgrounds in a cr-model context which naturally 
appear when one considers a quantum field theory of D-particles flT7| , |2"6fl . Unfortunately, the 
precise form of the effective moduli space potential in that case is still not known in a closed 
form (see [ID] for the first few terms in a velocity expansion of the C-function). For instance, 
it is known [0,|3lJ that when a string stretches between two D-particles one obtains a linear 
(confining) potential, but the inclusion of an arbitrary number of strings complicates the 
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situation. Moreover, as discussed in section V, there is no unique way of associating the 
non-linear Schrodinger equation ( p.6[ ) to the effective Fokker-Planck equation for the proba- 
bility distribution (|4.19 ) which was derived by world-sheet renormalization group methods. 



In particular, the precise form of the effective potential in the linear part of the associated 



Schrodinger equation is irrelevant for yielding the Fokker-Planck equation ( 4.19|) . This latter 
equation is essentially that which governs the quantum mechanics on moduli space. As we 
now show, the present approach enables one to make some non-trivial (and generic) state- 
ments about the existence of bound state solutions in the multiple D-particle system under 
consideration, without precise knowledge of the effective moduli space potential Vm(|7, u',t). 
To this end we discuss a second class of stationary solutions of the non-linear Schrodinger 



equation (|5\6| ), which are valid for a generic (but appropriate) potential Vm- Such a case 
incorporates the complex dynamics that binds strings and D-particles in "fat D-brane" 
configurations. Let us make the ansatz 

9 2 (y,t) = -mriMy) - ut (6.7) 



In the Galilean invariant case ( |5.11| ) of relevance to us here, the function ( |6.7p also solves 



( p. 3D and the resulting phase equation for 9i is 

2^V^x + 4(r? 2 + m ) (V,A) 2 - ^ = ^ (6.8) 

Then the function 

V (y)=expf36 1 (y) (6.9) 

where 



^ = 1 - — 
rj 2 J 4a'l 

with T> the diffusion coefficient ( 4.13| ), satisfies the affiliated linear Schrodinger equation 



/3 = 21 + -^ =1 (6.10) 
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2m{j} 2 )' 



V„.y + Vm (p = -utp (6.11) 



V2+V5 

Eq. ( |6.11|) corresponds to an ordinary linear Schrodinger wave equation of a particle in the 
potential Vm, but with a shifted effective mass 

m* = (3m (6.12) 



The stationary solutions of the non-linear equation ( |5.6| ) in the case ( |5.11fl are therefore 
given by 
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*(y,t) = <p(y)^f^ exp [~ J- gA _ p2 log^(y) 2 - iwfj (6.13) 

The wavefunctions ( |6.13| ) are square-integrable, \1/ G L 2 (M), if ^(y) 1 ^ are. 

Depending on the form of the potential Vm, such solutions may correspond to bound 
states of the affiliated linear equation ( |6.11| ), with the same potential Vm but with a shifted 
effective mass ( p,12| ). It is desirable that the effective mass be positive, or otherwise the 
interaction potential changes its sign. This requires j3 > in ( |6.10| ). For a fixed string 
coupling constant, from ( [4.13| ) we see that the requirement of positivity of the effective 
bound state mass m* translates into a limiting velocity for the D-brane dynamics. In the 



abelian case, this limiting velocity is just the speed of light f36fl , as follows from the form 
of the abelian Born-Infeld action which coincides with the standard relativistic free particle 
action in this case. Here we find the bound on the velocities of the multiple D-particle 
system appropriate to the non-abelian Born-Infeld action and the formation of bound state 
composites. Alternatively, for a given order of magnitude of the diffusion coefficient (f4. 13|) 
set by the recoil velocities of the D-particles, we find that bound state solutions exist for 
string couplings stronger than a critical value g* given by 



fj. 



c 9 _ 

°Y, tT \ ui \ 2 ( 6 - 14 ) 



\ 2 

The solutions of the non-linear Schrodinger equation differ from the bound state solutions 
of the affiliated linear equation by a spatially dependent phase proportional to the diffusion 
coefficient T>. In the framework of the "fat brane" picture of P^| , such bound states may be 
thought of as corresponding to bound state solutions of a system of D-particles interacting 
via the exchange of strings. The existence of a critical coupling (|6.14|) is then physically 
appealing, because it implies the formation of bound state condensates of D-particles and 
strings for strong enough string interactions, in similar spirit to the conventional quantum 



field theoretic case. This is in agreement with the fact |35| that for weak string interactions 
between the constituent D-branes, the effective spacetime is commutative and no bound 
states can form. When the string interactions become strong, bound states form and render 
the effective target space geometry noncommutative. Note that in the present case, the 
critical coupling (|6.14|) depends non-trivially on the (kinetic) energy scale of the recoiling 
D-particles, indicating that a certain minimum amount of energy is required to form the 
bound state. This velocity dependence is also quite natural, since it means that no bound 
states can form if the D-particles move too quickly relative to one another. 

Let us note also that QS.13j ) is the analog of the Haag decomposition in the non-diffusive 
case, since it is obtained from a linear Schrodinger equation via a non-linear, quantum 
mechanical gauge transformation. This non-linearity implies that the bound state solutions 
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of the affiliated linear wave equation with shifted mass m* and potential Vm.(v,u) acquire 
a phase that depends on y through the Haag amplitude <f(y}- The transition from the 
non-linear equation to a related linear one is "mediated" in part by this ^-dependent phase. 
The effective bound state mass ( |6.12| ) is then the analog of the Bargmann superselection 
criterion in the present case, which we see is modified by a non-linear transformation of the 
original BPS mass m and also the quantum recoil diffusion. 



C. Gaussian Wave Solutions 

The non-linear Schrodinger equations admit additional non-trivial solutions of a Gaussian 
wave type. Such configurations acquire particular importance in the "fat D-particle" context, 
given that among them there are solitary waves, and one would expect their appearance in 
view of the role of D-branes as solitons in string theory. The Gaussian wave ansatz for the 
wave-functional \l/ is described by 

t) = ~ (Mt)y 2 + B(t)y + C(t)) (6.15) 

where we have restricted ourselves, for simplicity, to one spatial dimension y, given that 
the Gaussian cases are separable [pTfl . Here s, a, A, B, C are real- valued functions of t which 
are determined from the amplitude and phase equations ( |6.3|) by equating the different 
powers of y that appear. This leads, in general, to a system of coupled non-linear ordinary 



differential equations for these functions of t [31,37|. This set of equations can be reduced to 



two ordinary second order differential equations for a and s, and A, B, C may be expressed 
directly in terms of the solutions of these equations. 

For our purposes we note that one may explicitly construct Gaussian wave solutions in 
the free case Vm — 0. Among them are also Gaussian solitary wave solutions with time- 
independent width a(t) = a = const. For the free case Vm = 0, the amplitude and phase 
equations (|6.3|) imply that a necessary condition for the existence of such solutions is 



V2 + V5 = a'f s - 7 ^- = (6.16) 



which, on account of ( |4.13j ) and ( |4.15| ), would imply extreme fine tuning between the mag- 



nitude of the velocities of the recoiling D-particles and the string coupling g s . Thus it is 
unlikely to be met in a generic situation, so that the class of Gaussian solitary waves do 
not seem to describe the D-brane dynamics (see however the discussion below for a different 
interpretation). 
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Nonetheless, the Gaussian ansatz (|6.15 ) allows one to carry out some explicit calculations 
and thereby see directly the effect of the quantum recoil through the non-linearity of the 
wave equation. This also enables one to make some non-trivial consistency checks of the 
present formalism. As an illustration, consider the free case Vm = 0, for which the differential 
equation for the function s(t) is that of the classical motion of a free particle, s = 0. Galilean 
invariance may then be used to set s(t) = for all t, because a general solution s(t) = u t+so 
can be obtained from the s = one by a G(l) transformation. The remaining differential 



equations can be easily integrated by quadratures to give the solutions |37 

a(t) = 



A(t) 



B(t) 
C(t) 




(6.17) 



where a , f and C are appropriate integration constants and ft is defined in (|6.10|) . The 
complete solution for the wavefunction \1/ is then obtained by substituting (|6.17| ) into (|6.15 ) 
and (|6.1|) with s = 0. From this solution it is straightforward to compute expectation values 
of operators in terms of Gaussian integrals, and in particular arrive at the variance relation 



AyAp = 2g s 



\ 



1 + 



+ fo 



V 



(6.18) 



The uncertainty relation ( |6.18|) corresponds to that of a minimum uncertainty 
wavepacket, with effective Planck constant determined by the string coupling, as is conjec- 
tured in |nj. In fact, eq. ( 6.18|) has a remarkably similar form to the modified Heisenberg 
uncertainty relations that were derived in |I| directly from the string er-model genus ex- 
pansion, in that it contains a diffusion term proportional to the total kinetic energy of the 
D-particles and their open string excitations which implies that the accuracy with which 
one can measure position and momentum in a system of D-particles depends on the (recoil) 
energy that arises in the measurement process. Such an energy dependence is characteristic 
of string-modified uncertainties and quantum gravitational effects. At time it = we can 
expand the right-hand side of ( ETTB ) for slowly moving branes to get 



AyAp~g s (2 + f ( 



^]Ttr|*7f 



(6.19) 



32 



To leading orders in cr-model perturbation theory and for weakly coupled strings, the mo- 
mentum fluctuations Ap and recoil velocities U coincide up to the BPS mass of the D-branes 
pTjj] . Choosing sgn/ = ~ 1, it follows that (|6.19 ) has the usual form of the string-modified 
Heisenberg uncertainty relation: 

AyA P ^- 9s (2 + f^ + ^^(A P r (6.20) 

The solution to the equation for the extrema of Ay in (|6.20| ) as a function of Ap leads 
to the usual minimum measurable length, and hence to the standard expectations that it 
is not possible to probe distances smaller that the intrinsic string length or the spacetime 
Planck length using only D-particle probes (Note that for /o > there is no bound on the 
measurability of lengths in the spacetime). Notice that one may in fact use the complete form 
of the phase space uncertainty relations derived in [Kj, along with the modified Ehrenfest 
relations ( [5.12 ), as a necessary criterion to determine in general if a candidate wavefunction 
represents the true quantum states of the D-particles and strings. We shall not do so here, 
but rather carry on with the heuristic analysis and present one final possible description of 
the quantum dynamics of the system. 



D. Solitary Wave Solutions 



There exist other types of solitary waves which may occur for certain regions of the 
two-dimensional parameter space spanned by 1] 2 and 775. The construction of such solutions 
in the free case Vm = may be carried out as follows For the Galilean invariant case 
under consideration, if 9\ and 9 2 are solutions of the equations (|6.7|) and (|6.8|), then the 
Galilean boosted amplitude and phase 



9 1 (y,t) = 9 1 (y-ut) 

h(v,t) = 9 2 (y-ut,t) -mu-y 



mu ' 



(6.21) 



with u a constant velocity, are non- stationary solutions of (|6.3| ). Such transformations take 
the original Galilean invariant solution of the non-linear Schrodinger equation (|5.6f) ,( |5.11D ) 
with 774 = 0, 773 = —1, to solutions in the more general case 



- 277! (1 +Vs) +m = 
The resulting solution for the wave-functional \I/ then reads 



V(y,t) = ^0 exp 0i (y- ut) - 1 



Amr)\9\ (y — ut) + mu ■ y + 



mu ■ 



rj 3 t + cut 



(6.22) 



(6.23) 
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where \l/o is a constant. The associated probability density V = W^S is itself a solitary wave, 
moving with constant speed and without changing its shape: 

V(y,t) = |^ | 2 exp29 1 (y-ut) (6.24) 

In the particular case where 772 (772 + Vs) < 0, i.e. (3 < 0, one obtains square-integrable 
solitary wave solutions of the form ]3~T| 

toff 

4a' of -D 2 



*(jT, t) = ^0 cosh [k- (y-ut) 



x exp 



log cosh (k ■ (y — ut) \ — s -— k 2 t 



Notice that the condition (3 < is opposite to that for the existence of positive effective 
mass bound states in the affiliated linear Schrodinger equation (|6.11 ). This implies that 
such coherent solitonic states cease to exist for string couplings stronger than the critical 
value ( |6.14|) , for given recoil velocities of the D-particles. 

Such solitary waves may play an important role in the underlying dynamics of a sys- 
tem of D-particles. The fact that they arise in the "fat brane" picture indicates that the 
latter degree of freedom constitutes an effective coarse-grained description of macroscopic 
quantum coherent states that may characterize (under the specified conditions) a certain 
phase of the ground state of a system of D-particles interacting among themselves via the 
exchange of strings. From this point of view, the critical coupling situation ( J6.16 ) where 
g s — g*, corresponding to Gaussian solitary waves, may be thought of as a critical regime 
separating the bound state phase from the solitary wave coherent phase. It should be noted, 
however, that although the wavefunctions ( |6.13|) and ( |6.25|) are non-analytic around g s = g*, 
this change of solution should not be regarded as some phase transition in string coupling 
constant space. Rather, when the string interactions become weak enough to untighten the 
bound states of D-particles and strings, the quantum dynamics is described by minimum 
uncertainty wavepackets which produce the appropriate smearing to the quantum spacetime. 
Outside of this regime the D-particles behave as freely propagating solitary waves lending 
another interpretation to the notion of D-branes as string theoretic solitons. It should also be 
stressed that this heuristic picture is complementary to the picture discussed in |IU|| , where 
the microscopic structure of groups of recoiling D-particles under the exchange of strings 
has been studied in detail in a cr-model context. In the present picture one finds an effective 
quantum mechanics for a "coarsed grained" description of the multiple D-particle system, 
in which "lumps" of D-particles are described as single bodies obeying a non-linear diffusive 
Schrodinger dynamics which encodes their interactions. As mentioned before, this analysis 
is at present preliminary and must thought of as a heuristic picture for the D-brane dynam- 
ics, given that it neither incorporates supersymmetry nor the precise form of the effective 
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moduli space potential Vm- The above discussion has only centered around the free case as 
an illustration, which is correct only to a certain approximation, since the abelianization of 
the collective coordinates of the D-particles, on which the "free" particle picture is based, is 
not accurate beyond this mean field approximation. 



VII. CONCLUSIONS 

The above discussion completes the first analysis concerning the construction of a 
Schrodinger-type wave equation (here non-linear) for the wavefunctional of a collection of 
D-particles. This equation is the present proposal for the complicated quantum dynamics of 
Dirichlet branes, at least within the framework of worldsheet cr-model perturbation theory. 
Further analysis should concentrate on the applicability of such studies to interferometric 
devices, which might probe Planckian physics. The non-linear Schrodinger dynamics that we 
have described above to characterize generic systems of D-branes should also be compared to 
the linear modifications of Schrodinger wave equations in open quantum mechanical systems 
interacting with a stochastic type environment. Such types of interactions are customarily 
assumed in many approaches to quantum gravity. 

There are, in general, arguments against non-linear Schrodinger wave equations, mostly 
stemming from the fact that non-linearities may lead to superluminal propagation, i.e. mo- 



tion which goes faster than the speed of light (see the discussion in and references 
therein). In the D-brane picture advocated above, there is the advantage of expressing the 
system in a "closed" but non-linear form, i.e. writing a self-consistent quantum equation 
which comprises only known degrees of freedom, including recoil. In the above analysis we 
did not use environment operators, which are usually unknown and assumed to be generi- 
cally of a stochastic form. In the present approach we have demonstrated the existence of 
quantum diffusion and of stochasticity [p~Q] , [T~5l| in the sense of Gaussian probability distribu- 



tions obtained explicitly as a result of summing over all worldsheet genera (i.e. coming from 
the quantum string theory). 

Moreover, as shown in ||10|| , the resulting off-shell target space action underlying the 
Schrodinger dynamics near a fixed point in moduli space (where the entire approach is 
valid) is the non-abelian Born-Infeld action, which implies a limiting velocity (equal to 
the speed of light) for the D-particle dynamics ||36|| . This feature prevents superluminal 



propagation. Thus the issues of superluminal propagation due to a non-linear Schrodinger 



wave equation, which were originally advocated in PBfl , cannot apply to the D-particle 
cases. These issues express such conceptual difficulties as the possibility of using long- 
ranged quantum correlations to send instantaneous signals, and of communicating among 
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different quantum mechanical worlds in a "multiple world" interpretation. As noted in p8 



since non-linear gauge transformations on the wavefunctional can be found which linearize 
the non-linear Schrodinger wave equation, such claims cannot be correct in general. The 
Fokker-Planck equation that we have derived for the probability density may correspond to 
a linear Schrodinger wave equation, but of an open system, i.e. including interactions with 
environmental degrees of freedom (in agreement with the analysis of jT7L[5S]]). However, the 
non- linear wave equation involving only functionals of and thereby closing the system, is 
preferred. We do not need to appeal here to the unknown environmental degrees of freedom 
to describe the recoil-induced diffusion of D-brane dynamics. 

Although the arguments for superluminal propagation are overcome by the existence of a 
Born-Infeld Lagrangian for the DO-brane dynamics, a more direct approach is desired since 
the Fokker-Planck diffusion equation implies that an external environment is present for the 
D-brane dynamics. This may be argued through the fact that the non-linear effects that we 
have considered cannot produce superluminal propagation because the present formalism 
is self-consistent only for weakly-coupled strings with g s — > (corresponding to heavy D- 
particles and hence very small velocities). Light string states are then obtained by ^-duality 
under which the resulting Schrodinger wave equation may indeed become linear. 

Notice also that the time-reversal violation of the effective wave equation that we have 
obtained is implied directly (through entropy production) by the fact that decoherence is 
present in the system of D-particles. However, if CPT is conserved in a standard string 
theory, then T violation should occur if CP symmetry breaking is present. In the Born- 
Infeld Lagrangian one may have induced CP violation in four dimensions, for example, by 
the introduction of topological instanton terms. In string phenomenology such terms arise 
quite naturally upon compactification of the target space. The present approach which 
allows decoherence (due to the energy dependence of the quantum uncertainties) therefore 
provides a natural explanation of the supression of CP symmetry in string theory. But 
generally in decohering quantum mechanical systems CPT is not conserved, so the above 
approach also differentiates between standard particles and recoiling D-particles. 

We close by mentioning that non-linear quantum dynamics may be desirable from another 
point of view discussed recently in ||39|| . The non-linear evolution of correlated systems (i.e. 
those involving the dynamics of many "particle" states) may be free from the contamination 
by negative probabilities which inflict the corresponding linear case. The D-brane systems 
we are considering here always involve many D-branes interacting with each other via the 
exchange of closed string states. In the present picture [[H]] single D-particles are only 
viewed as a limiting case of a multi-brane system where the constituent branes are separated 
at macroscopic distances compared to the string length \/~a' . The non-linear dynamics, 
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therefore, which we proposed in this article, may be essential for guaranteeing complete 
positivity of the system in the generalized sense discussed in |J9| . If one accepts the modern 
viewpoint of string- inspired phenomenology of the observable world [^(J , according to which 
the latter is represented as a D-brane interacting with other branes via the exchange of closed 
string bulk states only, then the above-mentioned issues may be of crucial importance for a 
consistent formulation. 
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APPENDIX 

In this Appendix we present a world-sheet formalism in support of the identification of 
the world-sheet partition function of a Liouville string with the probability density in moduli 
space rather than a wavefunctional, as happens in the case of critical strings. 

We commence our analysis by considering the correlation functions among vertex opera- 
tors in a generic Liouville theory, viewing the Liouville field as a local renormalizat ion-group 
scale on the world sheet [p2fl . Standard computations yield for an A-point correlation 



function among world-sheet integrated vertex operators V, = J d 2 zVi(z,z) : 

A N =< Vi,... V iN > M = T(-s)fi s < ( J d 2 z^e a,l> ) s V n . . .V lN >^ =0 (7.1) 

where the tilde denotes removal of the Liouville field zero mode, which has been path- 
integrated out in ( |7. 1| ) . The world-sheet scale fi is associated with cosmological constant 
terms on the world sheet, which are characteristic of the Liouville theory. The quantity s is 
the sum of the Liouville anomalous dimensions of the operators Vi 



N 



oii Q Q 



1=1 
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The T function can be regularized [f42| , |22"|| (for negative-integer values of its argument) by 
analytic coninuation to the complex-area plane using the the Saaschultz contour of Fig. [I]. 
Incidentally, this yields the possibility of an increase of the running central charge due to the 
induced oscillations of the dynamical world sheet area (related to the Liouville zero mode). 
This is associated with an oscillatory solution for the Liouville central charge near the fixed 
point. On the other hand, the bounce intepretation of the infrared fixed points of the flow, 
given in refs. [4~2,2I|, provides an alternative picture of the overall monotonic change at a 
global level in target space-time. 

To see technically why the above formalism leads to a breakdown in the interpretation of 
the correlator An as a target-space string amplitude, which in turn leads to the interpretation 
of the world-sheet partition function as a probability density rather than a wave-function in 
target space, one first expands the Liouville field in (normalized) eigenf unctions {<p n } of the 
Laplacian A on the world sheet 

(j)(z, Z) =J2 C n<Pn = C O 0O + J2 </>n 00 °C A^ (7.3) 

with A the world- sheet and 

A(j) n = -e n (j) n n = 0, 1, 2, . . . , e = (0 n , m ) = 5 nm (7.4) 

The result for the correlation functions (without the Liouville zero mode) appearing on the 
right-hand-side of eq. ( |7. 1| ) is, then 

A N oc / U n ^ dc n exp(—^- ^ e n c 2 n - ^ R n c n + 

E Oi0»(*)c»)( / d 2 ^e a ^^Y (7.5) 



with R n = J d 2 ^R^ 2 \^)(f) n . We can compute (|7.5|) if we analytically continue |4l| s to a 
positive integer s — > n 6 Z + . Denoting 

/(,,„) s •£ mtM (7 . 6) 

n,m ^0 €n 

one observes that, as a result of the lack of the zero mode, 

Af(x,y) = -4n6 (2 \x,y)-j (7.7) 

We may choose the gauge condition J d 2 £y/^<j) = 0. This determines the conformal properties 
of the function / as well as its 'renormalized' local limit 

f R (x, x) = lim x ^ y (f(x, y) + \nd 2 (x, y)) (7.8) 
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where d 2 (x,y) is the geodesic distance on the world sheet. Integrating over c n one obtains 
A n+N oc exp[^ a i a j f(z i , Zj) + 

Q 2 r r n , s n , s r , ^ Q 



1287T 5 



R(x)R(y)f(x, y)-J2 ^ / y/^R(x)f(x, z t )] (7.9) 



We now consider infinitesimal Weyl shifts of the world-sheet metric, j(x, y) — > 7(2;, y)(l — 
<j(x,?/)), with y denoting world-sheet coordinates. Under these, the correlator trans- 
forms as follows 



SA N oc \S^hia{zi) + — - / d 2 xJ^Ra(x) + 

l07T 7 v 

{Qs f d 2 xJ^a(x) + (s) 2 f d 2 x\f^a(x)f R (x,x) + 



Q s J J d 2 xd 2 y^R(x)a(y)G(x,y) - sY^®i f d 2 x^a(x)Q(x, z { ) - 

i 

^sY,aif R (zi,Zi) J d 2 x\fia(x)- 

Qs 1 f d 2 xd 2 y^{x)^y)R(x)f R (x ) x)a(y)})A N (7.10) 



16/T 

where the hat notation denotes transformed quantities, and the function Q(x,y) is defined 
as 

g{z,u) = f(z,u) - ±(f R {z,z) + f R (u,u)) (7.11) 

and transforms simply under Weyl shifts [22]. We observe from ( [7.101 ) that if the sum of 
the anomalous dimensions s / ('off-shell' effect of non-critical strings), then there are 
non-covariant terms in (|7.10|) , inversely proportional to the finite-size world-sheet area A. 
Thus the generic correlation function A^ does not have a well-defined limit as A — > 0. 

In our approach to string time we identify |22] the target time as t = 0o — — l°gA, 
where 0o is the world-sheet zero mode of the Liouville field. The normalization follows 
from a consequence of the canonical form of the kinetic term for the Liouville field <p in 
the Liouville a model | 43|J22 |. The opposite flow of the target time, as compared to that of 
the Liouville mode, is, on the other hand, a consequence of the 'bounce' picture ^2 , 22] for 
Liouville flow of Fig. 1. In view of this, the above-mentioned induced time (world-sheet scale 
A-) dependence of the correlation functions A^ implies the breakdown of their interpretation 
as well-defined S'-matrix elements, whenever there is a departure from criticality s ^ 0. 

In general, this is a feature of non-critical strings wherever the Liouville mode is viewed 
as a local renormalization-group scale of the world sheet |]22|] . In such a case, the central 
charge of the theory flows continuously with the world-sheet scale A, as a result of the 
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Zamolodchikov c-theorem [f44| . In contrast, the screening operators in conventional strings 
yield quantized values Due to the analytic continuation curve illustrated in Fig. [I], we 
observe that upon intepreting the Liouville field (f) as time t oc logA, the contour of Fig. 
1 represents evolution in both directions of time between fixed points of the renormalization 
group: Infrared fixed point — > Ultraviolet fixed point — > Infrared fixed point. 

When one integrates over the Saalschultz contour in fig. [I], the integration around the 
simple pole at A = yields an imaginary part [j|2[22|, associated with the instability of 
the Liouville vacuum. We note, on the other hand, that the integral around the dashed 
contour shown in Fig. [l], which does not encircle the pole at A = 0, is well defined. This 
can be intepreted as a well-defined ^-matrix element, which is not, however, factorisable 
into a product of S— and S^ — matrix elements, due to the t dependence acquired after the 
identification t = — log A 

Note that this formalism is similar to the Closed-Time-Path (CTP) formalism used 



in non-equilibrium quantum field theories [45 . Such formalisms are characterized by a 



'doubling of degreees of freedom' (c.f. the two directions of the time (Liouville scale) curve 
of figure [l], i n each of which one can define a set of dynamical fields in target space). As 
we discussed above, this prompts one to identify the corresponding Liouville correlators An 
with ^-matrix elements rather than S'-matrix elements in target space. Such elements act 
on the density matrices p = Tr^l^ >< rather than wave vectors > in target space of 
the string: p out = $p in (c.f. the analogy with the S'-matrix, \out >= S\in >). 

This in turn implies that the world-sheet partition function Z x l of a Liouville string at 
a given world-sheet genus x ; which is connected to the generating functional of the Liouville 
correlators An, when defined over the closed Liouville (time) path (CTP) of figure [I], can be 
associated with the probability density (diagonal element of a density matrix) rather than 
the wavefunction in the space of couplings. Indeed, one has 

Zx,lW\ = J^/^ZxA^g 1 } (7.12) 

where {g 1 } denotes the set of couplings of the (non-conformal) deformations, 0o ~ l n ^4 is 
the Liouville zero mode, and A is the world-sheet area (renormalization-group scale). If one 
naively interprets Z Xt i\<j>0t 9 1 ] as a wavefunctional in moduli space {g 1 }, ^[(j)o, g 1 ], then, in 
view of the double contour of figure |], over which Z x ^ is defined, one encounters at each 
slice of constant O a product of ^[00; fl ,/ ] 1 ^[ < / , 0j <? J ]> the complex conjugate wavefunctional 
corresponding to the second branch of the contour of opposite sense to the branch defining 
ty[4>o, g 1 }- This is analogous to the doubling of degrees of freedom in conventional thermal 
field theories [f4"5| ]. Such products represent clearly probability densities V[t, g 1 ] in moduli 
space of the non-critical strings upon the identification of the Liouville zero mode 0o with 
the target time t [2^]. 
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In the above spirit, one may then consider the (formal) summation over world-sheet 
topologies Xi an d identify the summed-up world-sheet partition function J2 X Zx,l[4>0i 9 1 } 
with the associated probability density in moduli space. In the case of D-particles, dis- 
cussed in this work, the moduli space coincides with the configuration space (collective) 
coordinates of the D-particle soliton, and hence the corresponding probability density is 
associated with the position of the D-particle in target space. We stress once again that the 
above conclusion is based on the crucial assumption of the definition of the Liouville-string 
world-sheet partition function over the closed-time-path of figure [l} As we demonstrate in 
the main text, the specific D-brane example provide us with highly non-trivial consistency 
checks of this approach. 

Before closing this appendix we would like to give an explicit demonstration of the 
above ideas for the specific (simplified) case of recoiling (Abelian) D-particles. We shall 
demonstrate below that, upon considering the non-critical cr-model of a recoiling D-particle 
at a fixed world-sheet (Liouville) scale c/>o = InA, and identifying the Liouville mode with 
the target time, the Euclideanized world-sheet partition function can describe a probability 
density in moduli (collective coordinate) space. 

To this end, let us first consider the pertinent a model partition function for a D-particle, 
at tree level and in a Minkowskian world-sheet E formalism: 

where e~ 2 ~ InA 2 = \\iA (c.f. ( |3.11|) ), on account of the logarithmic algebra []23|. In our 
approach e -2 is identified with the target time. This is why in ( [7.1 3|) we have not path- 
integrated over X°, but we consider an integral only over the spatial collective coordinates 
X l ,i = 1, . . .9 of the D-particle. The combination of cr-model couplings egf + gf - may 
be identified with the generalized (Abelian) position eY % of the recoiling D-particle ( |3.6| ). 
Notice that, since here we have already identified the time with the scale e~ 2 > 0, the step 
function in the recoil deformations of the cr-model (|3.7| ) acquires trivial meaning. We shall 
come back to a discussion on how one can incorporate a world-sheet dependence in the time 
coordinate later on. 



Suppose now that, following the spirit of critical strings [13 1, one identifies the 
Minkowskian world-sheet partition function ( |7.13| ) with a wavefunctional ty[Y t ,<f) Q = t]. 
The probability density in Y l space, V[Y\t] = ^[Y\ t}^*[Y\ t], reads in this case: 

Z x=Q , L [Y\t]\ 2 = J DX l J DX^e-^/ E ^^+^ 

(| DXie^^^^-^-^/eE^W 9 ^-^) ®((y ^e-^/^Wwj ; (7.14) 

where X± = X % ± X' 1 . Upon passing to a Euclidean world-sheet formalism, and taking 
into account that the Yi independent factor can be absorbed in appropriate normalization 
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of the cr-model correlators, one then proves our statement that cr-model partition functions 
in non-critical strings can be identified with moduli space probability densities. 

Notice that similar conclusions can be reached even in the case where the time X° is 
included in the analysis as a full fledged world-sheet field and is only eventually identified 
with the Liouville mode. In such a case, by considering the probability density as above, 
one is confronted with path integration over X° = X ± X'° cr-model fields, which also 
appear in the arguments of the step function operators Q e (X±) in the recoil deformations 
( p.7|) , that are non trivial in this case. However, upon Liouville dressing and the requirement 
that the Liouville mode be identified with the target time, one is forced to restrict oneself 
on the hypersurface X_ = in the corresponding path integral / DX+DX°_ (. . .). As a 
consequence, one is then left with a world-sheet partition function integrated only over the 
Liouville mode X + = 20 (c.f. Z in (|7.12j )), and hence the identification of a Liouville string 
partition function with a probability density in moduli space is still valid, upon passing onto 
a Euclideanized world-sheet formalism. It can also be seen, in a straightforward manner, 



that summing upon higher world-sheet topologies, as in [10|, will not change this conclusion. 
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FIG. 1. Contour of integration in the analytically-continued (regularized) version of r(— s) for 
s £ Z + . The quantity A denotes the (complex) world-sheet area. This is known in the literature 
as the Saalschutz contour, and has been used in conventional quantum field theory to relate di- 
mensional regularization to the Bogoliubov-Parasiuk-Hepp-Zimmermann renormalization method. 
Upon the interpetation of the Liouville field with target time, this curve resembles closed-time-paths 
in non-equilibrium field theories. 
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